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1 Introduction 

The strict positive realness (SPR) of transfer functions is an important performance specification, 
and plays a critical role in various fields such as absolute stability/hyperstability theory [15, 20], 
passivity analysis [11], quadratic optimal control [3] and adaptive system theory [16]. In recent years, 
stimulated by the parametrization method in robust stability analysis [1, 4, 6], the study of robust 
strictly positive real systems has received much attention, and great progress has beem made [2, 5] [7]- 
[10] [12]-[14] [17]-[19] [21]-[35]. However, most results belong to the category of robust SPR analysis. 
Valuable results in robust SPR synthesis are rare. The following fundamental problem is still open 
[14, 18, 22, 23] [25]-[28] [30] [31]-[34]: 

Suppose a(s) and 6(s) are two n-th order Hurwitz polynomials, does there exist, and how to find a 
(fixed) polynomial c(s) .such that c{.s)/a{s) and c(s)/6(s) are both SPR? 

By the definition of SPR, it is easy to know that the Hurwitz stability of the convex combination 
of a{s) and b{s) is necessary for the existence of polynomial c(s) such that c{s)/a{s) and c{s)/b{s) are 
both SPR. In [13, 14, 19], it was proved that, if a(s) and b{s) have the same even (or odd) parts, such a 
polynomial c(s) always exists; In [2, 13, 14, 17, 22, 23, 25, 27, 30, 31], it was proved that, if n < 4 and 
a{s),b{s) e K {K is a stable interval polynomial set), such a polynomial c(s) always exists; Recent 
results show that [22, 23, 25, 27, 28, 30] [32]-[34], if n < 6 and a{s) and b{s) are the two endpoints of 
the convex combination of stable polynomials, such a polynomial c(s) always exists. Some sufficient 
condition for robust SPR synthesis are presented in [2, 5, 10, 17, 22, 23, 27], especially, the design 
method in [22, 23] is numerically efficient for high-order polynomial segments or interval polynomials, 
and the derived conditions are necessary and sufficient for low-order polynomial segments or interval 
polynomials. 

This paper presents a constructive proof of the following statement: for any two n-th order poly- 
nomials a{s) and b{s), the Hurwitz stability of their convex combination is necessary and sufficient for 
the existence of a polynomial c(s) such that c{s)/a{s) and c{s)/b{s) are both SPR. This also shows 
that the conditions given in [22, 23] are also necessary and sufhcient and the above open problem 
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admits a positive answer. Some previously obtained SPR synthesis results for low-order polynomial 
segments [22, 23, 25, 27, 28, 30] [32]- [34] become special cases of our main result in this paper. 



2 Main Results 



In this paper, P" stands for the set of n-th order polynomials of s with real coefficients, R stands for 
the field of real numbers, d{p) stands for the order of polynomial p(-), and H"- c P" stands for the 
set of n-th order Hurwitz stable polynomials with real coefficients. 

In the following definition, p(-) € P™, g(-) € P", /(s) = p{s)/q{s) is a rational function. 

Definition 1 [29] /(s) is said to be strictly positive real(SPR), if 



(ii) f{s) is analytic in Rc[s] > 0, (namely, q{-) £ ); 

(iii) Re[/0-c^)] > 0, Vco e R. 

If /(s) = p{s)/q{s) is proper, it is easy to get the following property: 

Property 1 [9] If /(s) = p{s)/q{s) is a proper rational function, q{s) G if", and Vu G 
P,Re[/(iw)] > 0, then p{s) e ff" U ff"-^. 

The following theorem is the main result of this paper: 

Theorem 1 Suppose a{s) = s" + ais"-i + • • • + a„ e ff", b{s) = s" + 6is"-i + ••• + &„ e i?", 
the necessary and sufficient condition for the existence of a polynomial c(s) such that c{s)/a{s) and 
c(s) /6(s) are both Strict Positive Real, is 



The statement is obviously true for the cases when n = 1 or n = 2. We will prove it for the case 
when n > 3. 

Since SPR transfer functions enjoy convexity property, by Property 1, we get the necessary part 
of the theorem. 

To prove sufficiency, we must first introduce some lemmas. 

Lemma 1 Suppose a{s) = s" + ais"~^ + • • • + a„ G iJ", then, for every k e {1, 2, • • • , n — 2}, 
the following quadratic curve is an ellipse in the first quadrant (i.e., > 0, i = 1, 2, • • • , n — 1) of the 
P"~-^ space {xi,X2, ■ ■ • ,a;„_i) ^: 



(i) d{p)^d{q); 



A6(s) + (1-A)a(s) efl'",Ae 



[0,1]. 



{4+1 ~ 4cfcCfc+2 = 0, 
Q = 0, 
le {l,2,---,n},l^k,k+l,k + 2, 



^When n = 3, the ellipse equation is (sec [27, 28] for details): 



(a2Xi — aiX2 - a-j)'^ - 4(ai - xi){asX2) = 0. 



When n = 4, the two ellipse equations are (see [25] [30]- [32] for details): 



{(12X1 + X3 — aiX2 — as)^ — 4(ai — xi){a3X2 — 02x3 — 040:1) = 0, 
040:3 = 0, 

(030:2 — 020:3 — 040:1)^ — 4(020:1 -|- 0:3 — 010:2 — 03)040:3 = 0, 
01 — xi = 0. 



When n = 5, the three ellipse equations are (see [33] for details); 
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and this ellipse is tangent with the line 

r ci = 0, 

\ le{l,2,---,n},l^k + l,k + 2, 

and the line 



ci=0, 

I G {l,2,---,n}Jj^k,k + l, 

n 

respectively, where c; := Yl {—^y~^''0'jX2i-j-i, / = 1, 2, ■ • • , n, ao = 1, a;o = 1, flj = if i < or i > n, 

3=0 

and a;^ = if i < or i > n — 1. 

Proof Since a{s) is Hurwitz stable, by using mathematical induction, Lemma 1 is proved by a 

direct calculation. 

For notational simplicity, for a(s) = s" + ais"-iH ha„ e iJ", b{s) = s" + 6is"-1h h6„ € i?", 

Vfc e {1,2, •• -,71-2}, denote 

ci = 0,1 e {1,2, - ■■ ,n},l k,k + l,k + 2}, 

and 

^efe •= {(2^i.a=2,---,a;„_i)| a!^^^ - 44dfe+2 < 0, 

(i( = 0,; e {1,2, •••,«},/ ^ l,A; + 2}, 

n n 

where C; := J2i-'^y^^"'3^2i-j-i,di := J2i-'^y~^^bj^2i-j-i,l = 1,2, - ■ ■ ,n, ao = l,bo = l,Xo = 1, 

3=0 3=0 

a, = and 6^ = if i < or i > n, and Xi=Oifj<Oori>n — 1. 

In what follows, {A, B) stands for the set of points in the line segment connecting the point A and 
the point B in the space (xi, X2, • • • , Xn-i), not including the endpoints A and B. Denote 

■■= {(a;i,a;2,---,x„_i)| {Xi,X2,- ■ ■ ,Xn-l) G Ur=l'i<i<n-2(^i' ^j)' 
V^iGn«„iG {1,2, 2}} 

and 

fl^ := {(a;i,a;2,-- •,a;„-i)| {xi,x2,-- ■,Xn-i) G Ur=fi<i<n-2(^»' -^j)' 

VBiGf^^,iG{l,2,...,n-2}}. 

Lemma 2 Suppose a(s) = s" + uis""! + • • • + a„ G ff", 6(s) = .s" + &is"-i + • • • + 6„ G i?", 
if ^ 0, take (xi, X2, • • • , x„_i) G 1^" n 0^ and let c(.s) := .s"-i + (xi - e)s"-2 + x2s"-3 + 

• • • + a;„_2S + (x„_i + e) (e is a sufficiently small positive number), then for — — and t-tt) we have 

a(s) o(s) 

e i?, Re[4^] > and Re[^$^] > 0. 

Proof Suppose (xi, X2, • • • , a;„_i) G n O'', let c(s) := s""^ + (xi - e)s"-2 + X2s"-^ H h 

Xn-2S + {xn-1 + e), £ > 0, £ Sufficiently small. 

When n = 6, the four eUipse equations are (see [34] for details): 

{a2Xi + X3 — aiX2 — 0,3)^ — 4(ai - xi)(as + 03X2 + aiX4 — X5 — a2X3 — a4Xi) = 0, 



{ 



aexi + 04x3 + 020:5 — 0314 — 0512 = 0, 050:4 — 043:5 — 0613 = 0, 06a;5 = 0, 

- 03x2 + 01x4 — X5 — 02x3 — 04x1)^ - A(a2Xi + X3 — a\X2 — a3)(a6a;i + 04x3 + 02X5 — 0314 - 05x2) = 0, 
x\ = 0, 05x4 — 042:5 — aeX3 = 0, 06^5 = 0, 
f {aQX\ 

\ a\ — x\ = 0, a2X\ + X3 — aiX2 — 03 = 0, 062:5 = 0, 

(05x4 — 04x5 — 062:3)^ — 4(o6a;i + 042:3 + 022:5 — 03x4 — 05x2)062:5 = 0, 

01 — XI = 0, 02X1 + X3 — 01X2 — 03 = 0, 05 + 03X2 + 01X4 — X5 — 02x3 — 04x1 = 0. 



(06X1 + 04X3 + 02X5 — 03X4 — 05X2)^^ — 4(05 + 03X2 + 01X4 — X5 — 02X3 — 04Xl)(05X4 — 04X5 — O6X3) = 



{ 
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Vw e R, consider 

a{juj) 

= I ,\i2 [cit^'^"-'^ + C2u;^("-^) + • • • + c„_ia;2 ^ 
I a{juj) \^ 

n 

where q := ^ (— ajX2i-j-i, ? = 1, 2, • • • , n, ao = 1, = 1, aj = if i < or i > n, and = if 
i<Oori>n — 1, and c{u)^) is a real polynomial with order not greater than 2(n — 2). 

In order to prove that Vw € R,'Re[ '^^j^^\ ] > 0, let f = w^, we only need to prove that, for any 
e > 0,s sufficiently small, the following polynomial fi{t) satisfies 

fl{t) := Cit"-1 + C2t"-2 + . • • + Cn-lt + c„ 

+£(r-i - c{t)) > 0, e [0, +oo). 

Since (xi, a;2, ■ • • , a;„_i) e il", by the definition of O", it is easy to know that 

gi{t) := cit"-^ + 02^-2 + • ■ • + c„_it + c„ > 0, Vi G (0, +oo). 

Moreover, we obviously have /i(0) > 0, and for any £ > 0, when f is a sufficiently large or sufficiently 
small positive number, we have fi{t) > 0, namely, there exist < ti < t2 such that, for all £ > 0, 

t e [0,ti] U [t2,+oo), we have fi{t) > 0. 

Denote 

Ml = inf gi{t), 

t&[tut2] 

Ni= sup \e-^-c{t))\. 

Ml 

Then Mi > and Ni > 0. Choosing < £ < — ^, by a direct calculation, we have 

fl{t) := Cir-l + 02^-2 + • • • + Cn-lt + C„ 

+£(i"-i - c{t)) > 0, ytG [0, +00). 

Namely, 

Va;ei?,Re[4^] > 0. 

Similarly, since (xi, 2:2, • • • , .t,j_i) £ fl^, there exist < < f4 such that, for all £ > 0, t G 
[Ojts] U [t4,+oo), we have /2(f) > 0, where 

/2 W := (ilt"^' + ^2*"-' + • • • + rfn-lt + dn 
n 

di :=J2{-'^y^'bjX2i-j-i, Z = 1,2, 

where bo = 1, xq = 1, bi = ii i < ov i > n, and Xj = if « < or i > n — 1, and d(a;^) is a real 
polynomial with order not greater than 2(n — 2) which is determined by the following equation: 

-£(ju;)"-^ + £ (-£) ,,2(„-l) , 
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Denote 

g2{t) := dxf-^ + (^2^-2 + • • • + dn-lt + dn, 
M2= inf 02 (t), 

tG[*3,t4] 

N2= sup \e-^-d{t))\. 

teltsM] 
Mo 

Then M2 > and N2 > 0. Choosing < £ < — , we have 

7V2 

Va;Gi?,Re[^^] >0. 

Ml M2 

Thus, by choosing < £ < min{— — , -rp}, Lemma 2 is proved. 

A"! N2 

Lemma 3 Suppose a{s) = s" + ais"-^ + • • • + a„ G iJ", b{s) = s" + 61 s"-^ + ■ • • + 6„ e if", if 
A6(s) + (1 - A)a(s) G i?", A e [0, 1], then n 

Proof If VA e [0, l],aA(s) := Xb{s) + (1 - A)a(s) G iJ", by Lemma 1, for any A G [0,1], 
^efc ~ 1' 2, • • • , n — 2, are n — 2 elhpses in the first quadrant of the i?""^ space {xi,X2, • • • , Xn-i)- 

VA G [0, 1], denote 

fl°'>- := {(xi,a;2,---,a;„_i)| (xi, 0:2, • • • , a;„_i) G [j7=ii<j<n-2i^Xi,Axj), 

V^A. G^^:^^e{l,2,•••,n-2}} 

Apparently, when A changes continuously from to 1, fl°'^ will change continuously from to 
fi'', and will change continuously from to f^^^, fc = 1, 2, • • • , n — 2. 

Now assume n fi'' = </>, by the definitions of and and Lemma 1, 3ui > 0, U2 > 0, ui ^ 
ai,ui 7^ 61, and 3k G {1,2, •• -,72 — 2}, such that the following hyperplane L in the space 

{xi,X2, ■ ■ ■ ,Xn-l) 

Xi X2 Xn-1 ^ 

Li : 1 h • • • H — i 

Ml U2 Un-1 

separates and O'', meanwhile, L is tangent with flg^, f2g2) ' ' " > ^e(n-2) ^'^^ ^efe simultaneously (or 
tangent with Og^, r2g2, • • • , ^e(n-2) ^'^'^ ^"fc simultaneously). 

Without loss of generality, suppose that L is tangent with 0^^, fJgj) ' ' " > ^e(n-2) ^'^'^ ^efe ^^^Ita- 
neously. 

In what follows, the notation [x] stands for the largest integer that is smaller than or equal to the 
real number x, and {y)^ stands for the remainder of the nonnegative integer y divided by the positive 

integer z ^. 

Since L is tangent with fi"^, rjg2, • • • , ^e(n-2) ^^'^ ^efc simultaneously, note that a{s) is Hurwitz sta- 
ble and ui > 0,ui ^ ai,U2 > 0, using mathematical induction, by a lengthy calculation, we know that 
the necessary and sufficient condition for L being tangent with Clgi,Q,^2^ ' ' " >^e(n-2) simultaneously 

^For example, [1.5] = 1, [0.5] = 0, [-1.5] = -2, and (0)2 = 0, (1)2 = 1, (11)3 = 2. 
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is 4 

% -\- 1 r ^ 1 r ^ 1 

±{-li-^\,ut^<?~^~^^ =0 (1) 

and 

«, = (-l)^ 2 ^4^>^U22 ^ i = 3,4,...,n-l, (2) 

where ao = 1. 

Since Uj = (—1) 2 ^^^z ^ j = 3^ 4^ . . . ^ n— 1, L is tangent with , by a direct calculation, 

we have 

i + I r*^] r^i 

j2i-li^h,u^^+'K[^^~^^=0 (3) 

where &o = 1- 

From (1) and (3), we obviously have VA e [0, 1], 

i + 1 r^i r*i 

±i~li^^a,.ut''^u?~ ^ =0 (4) 

i=0 

where a\i := ai + \{bi — a^), ag = 1, 5o = 1, i = 0, 1, 2, ■ • • , n. (4) and (2) show that L is also tangent 
with ri''^(VA G [0, 1]), but L separates fl"-- and ^I'^r, and when A changes continuously from to 1, Q"- 
will change continuously from fi^j^ to fi^^, which is obviously impossible. This completes the proof. 

Lemma 4 Suppose a{s) = s" + ais"-iH hfln e H'^,b{s) = s" + 6is"-1h \-bn G H'^,c{s) = 

s"-i + a;is"-2 + • • • + x„-i, if Vw G ii, Ref^^l > and Re[^$^] > 0, take 

a{ju)) b{ju}) 

c (s) := c(s) + 5 ■ h{s), 6 > 0, S sufficiently small, 

c (5) c ( 

(where h(s) is an arbitrarily given monic n-th order polynomial), then — -— and — -r^ are both 

a{s) b{s) 

strictly positive real. 

Proof Obviously, 9(c) — d{a) = n, namely, c(s) and a{s) have the same order n. Since a(s) G .ff", 

there exists wi > such that, for all I w |> wi, we have Re( ) > 0. 
a{3^) 

*When n = 3, we have (see [27, 28] for details): 

U1U2 — aiU2 — a2Ui + as = 0. 
When n = 4, we have (see [25] [30]-[32] for details): 

{uiu'2 — aiUj — a2UiU2 + a'iU2 + a4Ui = 0, 
U3 = —UlU2- 

When n = 5, we have (see [33] for details): 

— oi«| — a2UiU2 + a3U2 + a4Ui — as = 0, 

U3 = — tllU2, U4 = — ■"I- 

When n = 6, we have (see [34] for details): 

{MiUj — aiUj — a2Miu| + a3w| + a4UiU2 — 051*2 — oetii = 0, 
U3 = —U1U2, Ui = — «5 = 
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Denote 

M3= inf Re(4^), = sup | Re(^) |. 

|a;|<a;i a{jUj} a{jw) 
Ms 

Then M3 > and N3 > 0. Choosing < e < — — , it can be directly verified that 

Similarly, d(c) = d(b) = n, and there exists (J2 > such that, for all I a; |> uj2, we have Ref f^*^'^! ) > 

b(juj) 

0. 

Denote 

Then M4 > and A^4 > 0. Choosing < e < it can be directly verified that 

Re(^^) >0,Va;e J?. 

Thus, by choosing < £ < min{^^, -^}, Lemma 4 is proved. 

The sufficiency of Theorem 1 is now proved by simply combining Lemmas 1-4. 

Remark 1 From the proof of Theorem 1, we can see that this paper not only proves the existence, 
but also provides a design method. In fact, based on the main idea of this paper, wc have developed 
a geometric algorithm with order reduction for robust SPR synthesis which is very efficient for high 
order polynomial segments [26]. 

Remark 2 The method provided in this paper is constructive, and is insightful and helpful in 
solving more general robust SPR synthesis problems for polynomial polytopes, multilinear families, 
etc.. 

Remark 3 Our main results in this paper can also be extended to discrete-time case. In fact, by 

applying the bilinear transformation, wc can transform the unit circle into the left half plane. Hence, 
Theorem 1 can be generalized to discrete-time case. Moreover, in discrete-time case, the order of the 
polynomial obtained by our method is bounded by the order of the given polynomial segment [27, 29]. 

c( s] c{ S I CIS ) 

Remark 4 If ——^ and are both SPR, it is easy to know that VA G [0, 11, , , , ^ — ttt^ 

a{s) b{s) 'J L. " \a[s) + {1 - X)b{s) 

is also SPR. 

Remark 5 The stability of polynomial segment can be checked by many efficient methods, e.g., 
eigenvalue method, root locus method, value set method, etc. [1, 4, 6]. 



3 Conclusions 

We have constructively proved that, for any two n-th order polynomials a{s) and 6(s), the Hurwitz 
stability of their convex combination is necessary and sufficient for the existence of a polynomial c(s) 
such that c{s)/a{s) and c{s)/b(s) are both strictly positive real. By using similar method, wc can 
also constructively prove the existence of SPR synthesis for low order {n < 4) interval polynomials. 
Namely, when n < 4, the Hurwitz stability of the four Kharitonov vertex polynomials is necessary 
and sufficient for the existence of a fixed polynomial such that the ratio of this polynomial to any 
polynomial in the interval polynomial set is SPR invariant [22, 23, 25, 27, 30, 31]. The SPR synthesis 
problem for high order interval polynomials is currently under investigation. 
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